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Abstract. Chiral discrimination with high harmonic generation (cHHG method) 
has been introduced in the recent work by R. Cireasa et al (Nat. Phys. 11 , 654 
- 658, 2015). In its original implementation, the cHHG method works by detecting 
high harmonic emission from randomly oriented ensemble of chiral molecules driven 
by elliptically polarized field, as a function of cllipticity. Here we discuss future 
perspectives in the development of this novel method, the ways of increasing chiral 
diclrroism using tailored laser pulses, new detection schemes involving high harmonic 
phase measurements, and concentration-independent approaches. Using the example 
of the epoxypropane molecule CaHgO (also known as 1,2-propylene oxide), we show 
theoretically that application of two-color counter-rotating elliptically polarized laser 
fields yields an order of magnitude enhancement of chiral dichroism compared to single 
color elliptical fields. We also describe how one can introduce a new functionality to 
cHHG: concentration-independent measurement of the enatiomeric excess in a mixture 
of randomly oriented left-handed and right-handed molecules. Finally, for arbitrary 
configurations of laser Helds, we connect the observables of the cHHG method to the 
amplitude and phase of chiral response, providing a basis for reconstructing wide range 
of chiral dynamics from cHHG measurements, with femtosecond to sub-femtosecond 
temporal resolution. 


Submitted to: J. Phys. B: At. Mol. Opt. Phys. 



Opportunities for chiral discrimination using high harmonic generation in tailored laser fields2 

1. Introduction 

Ever since their discovery, chiral molecular systems puzzled and inspired researches. 

The mirror symmetry, characterizing two molecular enatiomers, is possibly the simplest 
broken symmetry to think of. The complexity and variety of interactions and 
occurrences in nature of different enatiomers, flowing from such a small difference in 
their structure, is amazing. Chiral nature of living matter poses both fundamental and 
practical questions, from the origins of homochirality in biomolecules [Tj to detection 
and manipulation [2] of chiral properties. The oldest method used to detect chirality 
in such media is optical rotation: the polarization plane of a linearly polarized light 
propagating through left and right chiral media rotates in opposite directions (see e.g. 
optical rotation measurement in epoxypropane molecule; we shall use this molecule to 
illustrate our concepts here 0 ), 

Microscopically, optical rotation is governed by the interplay of dipole transitions 
in a chiral medium caused by electric and magnetic fields of a light wave. While 
the electric transition dipoles are identical in the two enantiomers in their respective 
molecular frames, the magnetic dipoles are pointing in opposite directions leading to 
opposite dynamical electronic response in the two enantiomers. The same mechanism 
underlies other linear chiroptical techniques which use circularly polarized light to detect 
differences in absorption (circular dichroism), Raman scattering, or circular fluorescence 
[4] between the two enantiomers. It is useful to note that while chirality is a topological 
property of molecular structure, chiral response to light is caused by chiral electronic 
dynamics, which thus effectively ’measures’ this structural property. 

Weak magnetic effects lead to weak chiroptical response. For example, circular 
dichroism is often three to six orders of magnitude smaller than light absorbance at 
the same wavelength. Weakness of chiral response poses challenges for time-resolved 
measurements of chiral dynamics. 

One way of increasing chiral response is to employ techniques that do not rely on 
weak interaction with the magnetic field component of the electromagnetic wave, such 
as e.g. photoelectron spectroscopy [5], microwave detection [6J, or Coulomb explosion 
imaging [Tj. Two of these techniques, one-photon [HI El ID] and multiphoton puna mg 
photoelectron circular dichroism (PECD) are pertinent to the approaches to chiral 
discrimination involving high harmonic generation (cHHG method [Hj). In particular, 
one-photon ionization is the inverse of photorecombination, which is the key step in high 
harmonic generation. 

The PECD detection requires resolving the direction of the final electron 
momentum: the photoelectron angular distribution obtained from randomly oriented 
chiral molecules shows asymmetry with respect to the direction of light propagation. 

This asymmetry originates from the electron interaction with the chiral potential of the 
core [HE], i-e. is the dynamical consequence of the chiral structure. The effect is most 
prominent for low energy electrons [TE] . since interaction with the chiral core potential 
is maximized in this case. The chiral dichroism signal achievable with PECD, up to 
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some 10%, sets the ’golden standard’ for the cHHG method pH]. 

With the advent of ultrafast light sources, there has been an increasing interest in 
developing ultrafast chiroptical spectroscopic techniques to investigate chiral response 
in the time domain. In the condensed phase, several breakthroughs have been achieved 
PH US] using non-linear chiroptical spectroscopy [musing & In the gas phase, where 
isolated compounds can be studied, such measurements are a lot more challenging due to 
low density of samples, so that no ultrafast measurement at ~ 10 1 fsec time scale has yet 
been reported. High sensitivity of PECD may open a path towards such measurements, 
provided that ultrashort circularly polarized laser mmm or XUV pulses mm 
are used. 

Importantly, cHHG naturally offers sub-femtosecond temporal resolution and allows 
one to detect multielectron chiral dynamics on their natural attosecond timescale. The 
chiral response can be monitored with At ~ 0.1 fsec resolution m, orders of magnitude 
better than achieved so far (see e.g. [23]). Moreover, already the first experiments [14] 
demonstrated very high enantio-sensitivity of the cHHG approach, leading to about 2- 
3% chiral signal for nearly linearly polarized laser fields, with ellipticity e ~ 1%. Thus, 
cHHG may represent an exciting alternative to PECD for ultrafast studies in the gas 
phase. Here we discuss the opportunities for extending the first cHHG measurements 
[T4] to increase their sensitivity and versatility, striving to turn cHHG into a standard 
ultrafast chiroptical detection technique. 

High-harmonics generation (HHG) can be understood as a sequence of three steps 
m m- ionization in a strong laser held, laser-induced acceleration of the liberated 
electron, and its recombination with the parent ion, all within the same laser cycle. 

High harmonic generation spectroscopy can be viewed as a pump-probe technique. 
Ionization acts as a pump, simultaneously launching the continuum electron and the 
correlated multielectron dynamics in the cation. Recombination acts as a probe of such 
dynamics [261 EH [28]. The pump-probe delay is controlled by light oscillation with 
attosecond precision; the ionization and emission times can be accurately characterised 
[29, EH E2j- Effectively, the emitted light records a movie of the recombining system, 
with each harmonic representing a single frame mm- 

The ability to record multielectron dynamics on attosecond time-scale m m is 
one of the most exciting features of high harmonic spectroscopy. This multielectron 
dynamics can be visualised as the motion of a hole created in a molecule upon ionization. 

The shape, location, and momentum of the hole during recombination are mapped 
onto the properties of the emitted light - its amplitudes, phases, and polarizations 
[271 Ml S3, 135] . This link is crucial for cHHG. 

Indeed, one of the effective mechanisms of chiral sensitivity of HHG is directly linked 
to attosecond laser-driven dynamics of the hole between ionization and recombination. 
Equivalently, one can think about laser-driven multi-electron dynamics in the molecular 
ion, i.e. light-induced transitions between its electronic states. Chiral sensitivity arises 
due to the interference of two types of transitions between the electronic states of the 
cation m the electric dipole transitions driven by the minor component of the electric 


Opportunities for chiral discrimination using high harmonic generation in tailored laser fieldsA 

field, and the magnetic dipole transitions driven by the major component of the magnetic 
held. This inteplay of the magnetic dipole and electric dipole transitions occurs on the 
sub-cycle time-scale. The enantio-sensitive sub-cycle dynamics reflects the molecular 
structure, which dictates the strength of the strong-held driven chiral response. Note 
that this dynamics occurs with the nuclei still located at their positions in the neutral 
molecule: there is very little time for the nuclei to move during the fraction of the laser 
cycle between ionization and recombination. 

Optimal regime for the interference of the magentic dipole transition driven by 
the major component of the magnetic held and the electric dipole transition driven by 
the minor component of the electric held requires balance between the magnitudes of 
the respective interactions, \pH 0 \ ~ \edE 0 \. Here E 0 = H 0 are the amplitudes of the 
electric and magnetic helds of the laser pulse, /i is the matrix element of the magnetic 
dipole transition, d is the matrix element of the electric dipole transition, e is the 
cllipticity of the laser held, eE 0 is the amplitude of the minor component of the electric 
held. Since \p\ ~ 10~ 2 |d|, the balance is achieved for |e| ~ 1%, i.e. for nearly linear 
helds. Since /2 changes direction by 180° between the two enantiomers of the chiral 
molecule, constructive interference requires opposite values of helicity for right and left 
enantiomers, giving rise to chiral dichroism. Chiral dichroism at a few percent level for 
such almost linear helds observed in Ref. [T3] demonstrates high potential for chiral 
discrimination of this new technique. However, as with any new technique, its hrst 
implementation may not be optimal. 

The observed cHHG signal ra maximises for very low ellipticities e ~ 1%. On 
the one hand, the chiral response typically increases with increasing ellipticity, likely 
maximizing for circular rather than nearly linear helds. On the other hand, the harmonic 
signal drops quickly with increasing ellipticity. These conflicting requirements seem to 
preclude exploration of what could be a more favourable regime. Resolving this conflict 
should enhance sensitivity and hexibility of the cHHG technique. 

Tailoring driving helds to maximize the chiral response and optimize the cHHG 
is a clear direction forward. Specifically, one should aim to (i) maximise effects of the 
chiral hole dynamics in the molecular ion between ionization and recombination, without 
sacrificing high harmonic generation efficiency; (ii) explore application of other tools of 
high harmonic spectroscopy to chiral discrimination, complementing measurements of 
the harmonic spectra with measurements of the harmonic phase [2T, 28j [29], [30]. Below, 
we show how both aims can be achieved using tailored multicolor helds. 

The use of multicolor, shaped helds for cHHG should allow one to achieve two 
important goals. 

First, additional helds can couple ionic states populated during ionization to the 
ionic states not directly populated during the ionization step, but characterized by 
strong magnetic dipole transitions, for example due to the proximity of their respective 
orbitals to a chiral center. This would make the cHHG method truly general, applicable 
to many chiral molecules, even if the ionization step does not directly populate electronic 
states with strong magneto-dipole transitions. 
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Second, circular light has long been used to detect chiral response because such 
light has high chirality. However, the high harmonic yield maximises for linear fields 
and quickly drops with increasing e. Application of chiral multi-color fields resolves 
this problem. The family of counter-rotating two-color fields, introduced in [36], can 
combine circular or elliptic polarization with strong high harmonic response [ 22 ]. 

Here we consider the application of elliptic counter-rotating two-color laser fields, 

Ei = E 0 (ei cos(u>t)e x + sm(u>t)e y ), E 2 = E 0 (e 2 cos(2ut)'e x — sm(2ut)e y ), (1) 
to induce, enhance and manipulate the chiral HHG response from randomly oriented 
gas of chiral molecules. In Eq.(l) E 0 is the held strength, ei, 62 are the ellipticities of the 
fundamental held and its second harmonic, e x y are the unit vectors in the poarization 
(x,y) plane of the laser held. 

Below, we theoretically show that application of such helds should allow one (i) 
to improve the sensitivity of the cHHG method to rival the sensitivity of the PECD 
measurements, without relying on using the driving light with very small ellipticity, 

(ii) to increase the range of harmonics for which the chiral signal is strong. We 
also show how phase measurements of the chiral high harmonic signal allow one to 
determine the amount of specific enatiomers (enantiomeric excess) in a mixture, in a 
way that does not require one to know the medium concentration. Finally, we present 
the theoretical background for the reconstruction of time-dependent chiral response for 
arbitrary ellipticities and configurations of the driving laser helds. 

2. Chiral discrimination based on HHG spectral and phase measurements 

In this section, we apply the approach we have already used in [ljj to estimate the 
high harmonic signal vs ellipticity of the held given by Eq. ([!]). We set the ellipticities 
of both helds to be equal, ei = 62 = e. We will consider an epoxypropane molecule 
CH 3 CHCH 2 O (also known as 1,2-propylene oxide), which has been used in the hrst 
cHHG experiments mi- 

There are two key cHHG channels in strong held ionization of epoxypropane. These 
channels are associated with leaving the ion either in the ground electronic state X 
upon ionization, or in the hrst excited state electronic state A. In the former case, the 
returning electron can recombine with the ion in state X, in the latter case, with the 
ion in state A. Higher excited electronic states are well separated in energy from these 
two, and do not contribute to ionization in strong mid-infrared laser helds. 

Recombination of the continuum electron with the parent ion in these states gives 
raise to direct HHG channels, XX and AA. Here the hrst letter marks the ionization 
channel and the last marks the recombination channel. Laser-induced transitions 
between the states A" and A of the ion open additional cross-channels [35]: XA and 
AX. For example, strong held ionization can leave the ion in the state X, but laser- 
induced transitions can excite the ion to the state A. The returning electron now has to 
recombine with this state: the hole created in the molecule during ionization is different 
from the hole before recombination. This results in the XA channel in HHG. 
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Importantly, AA and AX are chiral sensitive since they involve the interference of 
the laser-driven electric dipole and magnetic dipole transitions [TJ]. Onr calculations 
ra show that ionization and recombination create strongly preferred direction for the A 
channel, while the ionization-recombination angular dependence for the channel X has 
no such preferred direction. Also, the cross-channel XA associated with the ionization 
to the channel A" followed by recombination to the channel A is suppressed relative 
to the channel AX. Along the direction of maximal ionization for the channel A, the 
strength of the channel XA is 3 x 10~ 2 compared to the channel AX. 

These considerations allow us to develop a very simple model of the overall process 

CU: 

d JN) = d xx (TV) + d AA (A) + d XA ( N) + d AA - (iV) ~ 

~ d xx (AO + d AA (A) + d^ (A, & max ) = 

= d AA (AT) [[1 + R AX \e~ iEAxT W + R XX ] (2) 

The symbol d/j denotes N-tli harmonic component of the laser-induced polarization 
associated with channel IJ , i.e. ionization that leaves the ion in the state / and 
recombination with the ion in the state J. The overline denotes orientation-averaged 
quantities. The complex-valued relative amplitude of the cross-channel, R A x, is given 
by orientation-averaged amplitude of the laser-induced transition from the state A to 
the state X. This amplitude takes into account all laser-driven transitions, including the 
transitions driven by the magnetic field component due to the strong magnetic dipole 
which couples these states. The cross-channel AX is maximized in the direction given 
by the solid angle Q max that optimizes ionization into A and recombination into A" in 
the molecular frame. The contribution of the cross-channel is still averaged over all 
angles in the plane orthogonal to this direction. 

In the last equality in Eq.Q, we have written down explicitly the relative phase 
evolution between the channels AA and AA, with t(N) being the time-delay between 
ionization and recombination and E A xt(N ) the difference in the energies between 
the ionic states A and A". Relevant molecular properties including transition matrix 
elements are described in detail in | 12 ], where we have performed calculations for the 
single-color elliptical field. Here we use the same input and strategy and only change 
the single-color elliptical field to the field given by Eq. (1). 

According to the calculations of the ionization and recombination amplitudes [IT], 
the relative contribution of the A A channel, Rxx- is close to unity. We have established 
PU that in the single color elliptical field, due to the weakness of the chiral sensitive 
cross-channels (AA, AA), these channels can only leave their fingerprint in the high 
harmonic spectrum near a specific harmonic order, for which the main channels AA and 
AA interfere destructively. The phase term E A xt(N ) in Eq.(j2]) controls this interference. 

Following the model in Eq. (J2]) , to estimate chiral dichroism in the two-color 
elliptically polarized counter-rotating fields, we need to know the ionization ij(A) and 
recombination t r (N ) times for each harmonic order N. This allows one to calculate the 
chiral cross-channel: the amplitude of populating the state A at the time t r (N ) after 
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the ion was created in the state A at the time tfiN). These times are the well-known 
solutions of the saddle point equations for high harmonic generation in such two-color 
fields [36] . Using Eq.Q we can now evaluate the normalized high harmonic signal vs 
cllipticity ei = £2 = e (see Eq.(l)). 

Fig. 1 (left panel) shows the calculated chiral dichroism as a function of e = e\ = e 2 , 
for Eo=0.0239 a.u. (1=2 x 10 13 W/cm 2 ) and the carrier frequency u= 0.0224 a.u. 

(A = 2035 nm). High values of chiral dichroism are observed in a wide range of 
harmonic numbers, not only in the region of the destructive interference of the XX 
and AA channels, which is around harmonic 43. Stronger chiral dichroism in this region 
(about 80 %) as compared to results obtained in weakly elliptical fields (about 3 %) [[2] 
is due to the order of magnitude higher cllipticity (e rv_/ 0.14). 

There are two factors leading to this increase for relatively small ellipticities. Firstly, 
in the range e « 1, chiral dichroism is proportional to cllipticity: increasing cllipticity 
14 times will respectively increase the chiral dichroism. Secondly, for the field used here, 
the amplitude of the transition between the states A and X is almost two times higher 
for e ~ 0.14 then for the single-color field of the same cllipticity. 

For higher ellipticities, the disparity between the laser-induced transitions between 
left and right enantiomers increases and the destructive interference is no longer the 
decisive mechanism in revealing chiral effects. For example, chiral dichroism at the level 
of about 20% is found for harmonics 35-51 and cllipticity e\ — €2 — 1. Such high values 
of chiral dichroism are on par with values achievable with PECD method. 

The specific value of chiral dichroism for large ellipticities depends on many 
additional factors. One such factor is the achiral background. Achiral contributions 
to the harmonic signal are associated [T4j with the ’diagonal’ HHG channels, involving 
the cation staying in the same electronic state between ionization and recombination. 

We found that importance of these channels in cHHG can be controlled via the Stark 
shift in the two-color field Eq.(l). Specifically, most chiral molecules are polar and have 
permanent dipoles in the cationic states. Laser interaction with the permanent dipole 
induces the linear Stark shift. For intense laser fields, these shifts easily approach the eV 
scale. Such shifts play important role in laser-induced dynamics in the cation between 
ionization and recombination. Therefore, they can be used to control cHHG by shaping 
the laser field in the plane of polarization. This is precisely what is achieved by varying 
ti = e 2 = e in Eq.(l). 

Indeed, for a given electronic state of the cation J with a permanent dipole d J , the 
additional phase accumulated in this state due to its interaction with the laser electric 
field is: 

tstark* f dtE 0 {ef x {t)di + fy{t)d J y ). (3) 

du 

Here f x>y (t ) describe the sub-cycle temporal structure of the electric field in the x,y 
polarization plane, d J x and dy are the corresponding components of the permanent dipole. 

The amplitude djj of the achiral contribution to the high harmonic field, associated with 
this state, acquires the corresponding additional phase factor, djj oc e^stark. In strong 
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mid-IR fields fstark ^ 1- Thus, the phase term exp {ifftstark) affect the outcome of 
coherent averaging of the harmonic response over all molecular orientations, in every 
channel. 

Changing e in Eq.(3) changes exp [i(j>g tark \, which affects coherent averaging over 
molecular orientations. Fig. 1 (right panel) demonstrates the linear Stark-based 
cllipticity control of the amplitude of the achiral channel AA in epoxypropane, using 
the same theoretical approach as in Ref. [[13] . Clearly, we can manipulate the amplitude 
of this channel within large dynamic range, approaching one order of magnitude in 
amplitude for the highest harmonics. The channel XX in epoxypropane shows similar 
(but not identical) behaviour. 
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Figure 1. HHG in epoxypropane in two-color counter-rotating elliptically 

polarized fields. Left panel: Estimated chiral dichroism Q(N,e) = 2[Ys(N,e) — 
Yr(N, c)]/[Ys(N, e) + Yr(N, e)], where Ys t R is the harmonic yield for left (S) and right 
(R) molecules. The color scale is in percentage. Right panel: The amplitude of AA 
channel, coherently averaged in the direction orthogonal to the direction of maximal 
ionization of the epoxypropane. Bottom panel: Relative phase (between the left and 
right enantiomer) in radians for the chiral sensitive channel AX in epoxypropane, 
without the contribution of the achiral channels. All calculations were done for electric 
field amplitude Eo=0.0239 a.u. (1.24 x 10 8 V/cm), fundamental frequency w=0.0224 
a. u. (A = 2035 nm), ei = 62 = e. 
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3. Chiral discrimination and chiral dynamic imaging based on HHG 
spectral and phase measurements 

High harmonic spectroscopy gives access to many observables, including harmonic 
phases PH 123 E3- Chiral sensitivity of the harmonic phase has not been explored 
so far. Fig. 1 (bottom panel) shows the relative phase between the chiral-sensitive 
transition amplitudes in the cation, for left and right enantiomers. This relative phase 
is substantial, up to one radian, well within the accuracy of the experimental phase 
measurements m- Thus, cHHG can (and should) also aim at using phase measurements 
to detect enantiomers. In general, achiral channels may hide this phase. However, here 
we are assisted by the suppression of achiral background shown in Fig.l (right panel), 
which indicates that the enantio-senitive harmonic phase could be observable for large 
ellipticities of the driving two-color field, where suppression of the achiral background 
is strong. 


3.1. Detection of the amount of enantiomers in a mixture. 


The combination of the harmonic phase measurement with the harmonic spectral 
measurement brings the new opportunity: concentration-independent detection of 
enantiomeric excess in a mixture M of left S and right R chiral molecules (see 
Fig.2(a)). Polarization of the medium is proportional to number density. Therefore, 
accurate determination of enantiomeric excess in chiroptical spectroscopies requires 
precise determination of number density, which can be very challenging. In the 
following we show how one can perform the measurement of enantiomeric excess 
without any knowledge about number densities of the mixture and of pure samples. 
Suppose we have a mixture of chiral molecules with an overall number density Nm 
and the relative fractions of the left-handed and right-handed molecules ms and m#. 
Our goal is to find a = ms/mR , which maps into the enantiomeric excess (e.e.) as 
e.e. = (1 — a)/(l + a) x 25%. 

The first step is to measure the harmonic yield from the unknown mixture, 
Ym{N,€). Ideally, this should be done for all e, changing it from e = 0 to e = 1 
and then reversing the hclicity (denoted as negative e). At the very least, one should 
measure the yield for, say, e = 1 and e = —1, i.e. perform two spectral measurements. 

Next, we use these two measurements to obtain the number density-independent 
quantity such as the elliptic dichroism, 


ED m (N,6 ) 


2 Ym(N , e) — Ym(N , —e) 

Yr,s,m(N, e) + Yr,s,m(N : —e) 


(4) 


or the ratio of the two harmonic intensities 

Y m (N, +e) 


Pm(N,€) = 


Y M (N,-e)' 


(5) 


We can then compare these measurements to the harmonic yields W,s(W, e) from 
the samples of pure S (left-handed) or R (right-handed) molecules, also made number 
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Figure 2. Cartoon illustrating importance of the harmonic phase in measuring 
enantiomeric excess using cHHG. (a) A mixture of 70% right-handed and 30% left- 
handed molecules. The goal is to find a = 0.3/0.7 (b) Schematic representation of the 
harmonic response from three gas cells with pure left-handed molecules, pure right- 
handed molecules, and the unknown mixture. The overall number density of the 
mixture is Nm, the relative fractions of the left-handed and right-handed molecules 
are ms and m R , Pr,s,m are the polarizations of the media, d Rt s,M are the harmonic 
dipoles, (c) In a mixture, the harmonic dipoles from the left-handed and right-handed 
molecules add coherently, therefore the relative phase between them is crucial for 
determining a = ms/m R and the enatiomeric excess. 


density-independent by measuring the corresponding elliptic dichroisms 

ED rs (N , e) = 2 Yr ' S ^' +6 | ~ Yr ' s [ N ’ ~ 6 \ ( 6 ) 

’ 1 ' Y RtS (N,+e) + Y R , s (N,-e) 1 J 

and/or finding the number density-independent ratios 

r 2 (n + e ) , 

U ' s{ ’ } Y R , s (N,-eY ( } 

For definitiveness, we assume that fd R > 1, i.e. for the R molecules the harmonic yield 
for e is higher then for —e. 

Suppose we find that for the pure R-sample and, say, e = ±1, the corresponding 
elliptic dichroism is ED R (N,e = 1) = 10%, and for our mixture we have measured 
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EDm(N,c = 1) = 5%. We also know that for the racemic mixture the dichroism is 
zero. Can we conclude that our mixture is exactly half-way between pure and racemic? 
Fig.2 shows why we cannot make such conclusion: the harmonic dipoles from S and R 
molecules add coherently in the mixture. Thus, the total signal also depends on the 
relative phase dsn of the harmonic emission between the left-handed and right-handed 
molecules - the angle between the two harmonic dipoles in Fig. 2 (c). This phase can 
be measured using the two-slit harmonic phase measurement technique developed in 
In this case, the measurement should use far-held interference of the harmonic 
light from two gas cells with enantiomerically pure S and R molecules. The shift of 
the intereference fringes relative to the case of two identical cells determines the relative 
phase <f>sR in a concentration-independent way, with relative concentrations determining 
fringe contrast but not fringe positions. 

Once we have measured fsRj we can use it together with the measured (3 ^ and Am 
to find a = mg/m r. Simple algebra following addition of the harmonic dipoles shows 
that 

„2 = a 2 + 2af3 R cos 4> S r +Pr 

M Pr® 2 + 2 OL^R COS (j)SR + 1 

Here we have also taken into account that | ds(N, —e)| = |d^(iV, +e)|. 

Varying a = mg /tur between a = 0 (pure R-sample) and a = oo (pure S- 
sample), we see that /3m is squeezed between /3r > 1 and /3s = 1 /(3r < 1 , i.e., 

As = 1 /An < $m < Afi- 

Next, we re-write the above equation as a quadratic equation for a and, taking into 
account that As = 1/An < Am < (3r , find its positive solution: 

-Ail cos <I>sr{ Am - 1) + V(3 2 r cos 2 4>sr(P 2 m - l) 2 + (A% ~ Am)(A_|Am ~ x ) , q x 
a Wh -1) • U 

where all quantities depend on the harmonic number, N, and the ellipticity e. 

Clearly, the main challenge in this approach stems from the experimental signal 
to noise ratio, which may result in using small quantities with large error bars. This 
challenge can be addressed by taking advantage of the huge redundancy of the two- 
dimensional cHHG measurement, performed for various harmonics N as a function of e. 
Indeed, every observable used in Eq. S depends on two parameters: N , e. Therefore, 
the application of Eq. (J9]) to every harmonic and every ellipticity should yield exactly 
the same value of a(N, e) = mi/m 2 . 

The above scheme includes a set of reference measurements performed on pure 
samples. Their outcome can be used to characterize multiple samples containing 
unknown mixtures of left and right molecules, under the same experimental conditions. 
The additional measurements, specific to the mixture, are the HHG intensity 
measurements yielding elliptical dichroism or (3 2 M (N 1 e) for the unknown sample. 

Crucially, identification of a using Eq. @i is concentration-independent: one does 
not need to know the overall number density of molecules in the pure R-sample, and 
one does not need to know the overall number density of molecules in the mixture. 





Opportunities for chiral discrimination using high harmonic generation in tailored laser fieldsl2 


Alternative scheme for detecting enantiomeric excess. The above measurement scheme 
for detecting enatiomeric excess is not unique. One can employ two additional phase 
measurements, specific for the mixture, measuring the relative phase between harmonics 
emitted by the left-handed molecules and the mixture and the relative phase between 
the harmonics from the right-handed molecules and the mixture. We now consider such 
a scheme. 

Again, we prepare 3 gas cells: with enantiomerically pure left molecules S, right R 
molecules, and the unknown mixture M. Using the two-slit technique, we now measure 
three relative phases: (i) the relative phase between the harmonics generated from the 
left-handed and the right-handed samples fsm (h) between the left-handed molecules 
and the mixture 4>ms, and (iii) between the right-handed molecules and the mixture, 
4>mr ■ We now define number density-independent signals by normalizing the harmonic 
fields at e 7 ^ 0 to their values at e = 0 , 

Ps,r(c) 


Ps,r {e = 0 ) 
Next, we note that 

P M (e) 


d sA e ) _ n / n 
d s A' = o) = DsMc) 

mid s (e) + m 2 d R (e ) 


( 10 ) 


= m 1 D s (e) +m 2 D R (e) (11) 


Pm{z = 0) mid 5 (e = 0) + m 2 d R (e = 0) 
since mi + m 2 = 1 and d s (e = 0) = d R (e = 0). Here D RR are the normalized values of 
harmonic fields as defined in Eqs. ( TofTT ). 

We now show how a = m\jm 2 can be extracted from the measurements described 
above. Writing D SR = |H Si u| exp(if>s,R), we can explicitly calculate the argument 
arg[DMD* R } = 4 >mr in two steps. First, 

D m D r = [miD s + m 2 D R ] D R = mxDsD^ + m 2 \D R \ 2 . (12) 

Second, for the quantity D S D* R in the above equation we write 

D s D r = \D s \\D R \mi [cos((Psr) +fsin(0 SjR )] +m 2 \D R \ 2 , (13) 

Using these two expressions, we finally obtain 


tan (4>mr) — 


m i\Ds\\D R \ sin(</>su) 


m i\Ds\\D R \ cos((/)sr) + m 2 \D R \ 2 ' 
Analogously, we obtain for 4>ms : 

tan (4>ms) = 


m 

2 \D : 

s 

\Dr\ 

sin (<j> SR ) 

m 2 \ 

l^sll 

Dr\ 

cos (4> sr ) + mi\D s \ 

2 


(14) 


(15) 


Introducing the ratio of the spectral amplitudes of the harmonic signals from left and 
right molecules P(N,c) = \D s (N,e)\/\D R (N,e)\ and j(N,e) = l |tn^Mg(v,e))p wllich 
can be obtained from the phase measurements, we can express a = m\jm 2 from Eqs. 

( PP| : 

,^ T \ cos(<j>sR(N, e))('y(N , e) - 1) + ^cos 2 (0 SH (A^, e))(l - i(N, e )) 2 + 4q(At, e)) 

“ (A, ' e) =-«- 16 
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Again, we note that every observable in Eq. (16) depends on N and e. Therefore, 


the application of Eq. (16) to every harmonic and every cllipticity should yield exactly 


the same value of a(N, e) = thus reducing the signal to noise ratio. 

We remind the reader that in this scheme, the independence of the measurement 
on concentration (number density) arises from normalizing the harmonic response to 
its value at e = 0, see Eqs.( TofTT ). From now on, all signals will be understood 
as normalized to their values at e = 0, thus explicitly removing their dependence on 
molecular concentration (number density). 

Note that in the second approach we perform redundant spectral measurements by 
measuring harmonic response for both S- and R-molecules while varying cllipticity in 
the full range, —1 < e < 1, increasing the redundancy of the first scheme by factor two. 


3.2. Time-resolving sub-femtosecond chiral response in full range of ellipticities 


As established in Ref.pHj, the main difference between the cHffG signals coming from 
the two enantiomers is due to the small magnetic-dipole response in the cation. Thus, 
the multi-electron wavepacket in the cation acquires small correction, which is linear in 
p = 1771 (a straightforward result of the first order perturbation theory): 

\*R, s (t)) = \Mt))±v\M*(t)). (17) 

ffere l'f'o(t)) is the time-dependent wavefunction describing the dynamics in the cation 
in the dipole approximation, ±//|5'& c h(t)) is the chiral correction of the time-dependent 
response, opposite for S and R molecules. This small correction to the multi-electron 
dynamics in the cation yields chiral correction to the harmonic field 

D r>s (N, e) = [D® q (N, e) ± pSD ch (N, e)]. (18) 

ffere we have used the number density-independent signals, which can be obtained by 
normalizing the harmonic fields at e ^ 0 to their values for e = 0, 


Ps,r(c) 


PsAe = 0 ) ds } n(e = 0 ) 


dsA e ) _ n / n 
= Us,R( e ) 


(19) 


ffere D^ S (N, e) is the normalized harmonic field in the dipole approximation. 

The first-order in p expression is general. Its only assumption is the linear response 
of the cation to the weak effect of the magnetic field. The dipole e) is not chiral 

sensitive (once averaged over randomly oriented ensemble), D^ y> (A, e) = D^\n, e) = 
D^ 0 \N,e). The chiral-sensitive component is p5D c ifN, e). The harmonic yield for left 
and right enantiomers is then: 

|AD ch (A0| 


Y R)S (N,e) = \DW(N)\‘ 


1 ± 2/i 


|ZT°)(A0| 


cos($ ch(N)) 


( 20 ) 


Using Eq. (J20 ), one obtains for the chiral dichroism Q(N,e ) = 2 (^ ie 

difference between the harmonic yields Ys t n(N,e) for left and right molecules divided 
by the signal from the racemic mixture): 


Q(iV,c) = 4f! cos($rf,(iV, e)). 


( 21 ) 













Opportunities for chiral discrimination using high harmonic generation in tailored laser fieldslA 


Here <f> c h is the relative phase between the chiral and achiral contributions to the emission 
amplitude. Thus, the standard measure of the chiral signal, when applied to high 
harmonic yield, is proportional to p,\SD c h(N, e)| cos (& c h(N, e)) - the chiral part of HHG 
response. 

An exciting goal is to separately reconstruct the amplitude /jj\8D^(N, e)| and the 
relative phase cos( < h c / l (A^, e)) of the chiral response. This can be done by augmenting 
the measurement of Q(N,e) with two more measurements: (i) of the relative harmonic 
phase 4>sr(N, e) between the S and R molecules, and (ii) of the harmonic amplitude in 
the racemic mixture |ZT°)(iV, e)|. 

Eq.(18) allows us to explicitly calculate tan(0s^): 

[■ D <°> - iiSD ch }[D^ + /j,5D ch \* « |ZA°)| 2 + ii5D* ch DW - p5D ch D^\ (22) 

where we have neglected the second order terms in //. Thus, 

\SD ch \ 


tan (4> S r) = 2/xsin($ ch )- 


(23) 


|ZT°)| ' 

Since the sign of becomes important, we specify that & c h = arg[AD* h ZT°)]. Adding 
tan(0g R ) 2 and Q 2 (N, e)/4 we obtain 


tan 2 (0 SR (AT,e)) + 


Q 2 (N,e) 


— W 2 (N , e), 


(24) 


where W(N,e) gives access to the relative amplitude of the chiral component: 

= < 25) 

The dependence of W(N, e) on | (N, e) \ can be further eliminated by performing the 

harmonic measurements in the racemic mixture (rm). Suppose this measurement yields 
the harmonic spectrum y rm (N,e ) (also normalized to the signal for e = 0). Then, the 
amplitude of the chiral component can be obtained as follows: 

W(N,e) y/y rm (N , e) 


pt\6D ch (N,e) | = 


(26) 


On the other hand, combining Eqs.(21), (25) we obtain the phase between the chiral 
and achiral components: 

Q(N,e) 


cos($ ch (N,e)) = 


(27) 


2W (N, e) 

Thus, both the amplitude and the relative phase of the chiral signal can be obtained 
from the combination of spectral and amplitude measurements. 

The result described by Eqs.( 25p6 27) is in fact more general then may appear 
at first glance: Eq.(18) may equally apply to chiral dynamics initiated in a neutral 
molecule and then probed with HHG. This scheme is sketched in Fig.3. Such pump- 
probe scheme could allow one to probe the evolution of chiral component p\5D ch (N, e) | in 
neutral species during structural changes on a longer, femtosecond time scale. This time 
scale is determined by the delay between the pump pulse, which excites chiral dynamics 
in the neutral molecule, and the probe pulse which generates the chiral-sensitive high 
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harmonic signal. Up to now, no pnmp-probe time-resolved measurement of ultrafast 
chirality has been done in the gas phase due to the absence of sufficiently sensitive 
technique. cHHG promises to open this direction. 


At 



Pump HHG probe 


Q(N,e, At)=2^ L J £ 0C fi (At) \s D ch Ov, e, At) | 

y l" i " y r 

Figure 3. A cartoon of pump-HHG probe measurement describing one of the possible 
schemes of detecting time-dependent chiral response. The chiral signal in HHG arises 
from the magnetic dipole. A pump pulse initiates nuclear motion in a chiral molecule. 
As a result of this motion, the nuclei rearrange themselves around the chiral center. 
This leads to changing g and therefore changes the high harmonic response induced 
by the HHG probe. Thus, HHG allows one to follow chiral response from sub fs to fs 
time scale. 


4. Conclusions 

Using the example of epoxypropane molecule, we have shown theoretically that cHHG 
method should strongly benefit from the application of the two-color, counter-rotating 
elliptically polarized fields. Application of such fields should lead to much higher values 
of circular dichroism compared to the original scheme ra. which relied on weakly 
elliptical fields. 

We predict that the combination of the harmonic amplitude and phase 
measurements should allow one to directly measure enantiomeric excess in a 
concentration-independent way and monitor the evolution of the chiral dynamics with 
sub-femtosecond time resolution. 

We have shown how cHHG method can be used to extract both the amplitude and 
phase of the chiral signal, for arbitrary configurations of laser fields used to initiate high 
harmonic process. Indeed, the analysis in Section 3 does not assume any specific field 
configuration. 
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We predict that Eqs.( 25p6p7 ) can also be used to extract the amplitude and the 
phase of the chiral signal from the HHG observables in the pump-probe set-up. In 
this setup chiral dynamics is first excited in a neutral molecule and then probed with 
chiral-sensitive high harmonic signal. From single shot to pump-probe, cHHG driven by 
tailored laser fields thus enables imaging of a wide range of chiral dynamics, with time 
resolution extended from sub-femtosecond to many tens of femtoseconds. 
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